arXiv:1501.05791vl [astro-ph.EP] 23 Jan 2015 


The Astrophysical Journal, 794:91 (12pp), 2014 October 10 

Preprint typeset using IATIjjX style emulateapj v. 5/2/11 


FORMING CHONDRULES IN IMPACT SPLASHES 
I. RADIATIVE COOLING MODEL 

Cornelis Petrus Dullemond, Sebastian Markus Stammler 
Institute for Theoretical Astrophysics, Heidelberg University, Albert-Ueberle-Strasse 2, 69120 Heidelberg, Germany 


Anders Johansen 

Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 43, 22100 Lund, Sweden 
The Astrophysical Journal, 794:91 (12pp), 2014 October 10 

ABSTRACT 

The formation of chondrules is one of the oldest unsolved mysteries in meteoritics and planet for¬ 
mation. Recently an old idea has been revived: the idea that chondrules form as a result of collisions 
between planetesimals in which the ejected molten material forms small droplets which solidify to be¬ 
come chondrules. Pre-melting of the planetesimals by radioactive decay of 26 A1 would help producing 
sprays of melt even at relatively low impact velocity. In this paper we study the radiative cooling of a 
ballistically expanding spherical cloud of chondrule droplets ejected from the impact site. We present 
results from a numerical radiative transfer models as well as analytic approximate solutions. We find 
that the temperature after the start of the expansion of the cloud remains constant for a time t coo i 
and then drops with time t approximately as T ~ To[(3/5)t/t coo i + 2/5] -5 / 3 for t > t coo \. The time at 
which this temperature drop starts f coo i depends via an analytical formula on the mass of the cloud, 
the expansion velocity and the size of the chondrule. During the early isothermal expansion phase 
the density is still so high that we expect the vapor of volatile elements to saturate so that no large 
volatile losses are expected. 

Subject headings: chondrules, radiative transfer 

1. INTRODUCTION 

The formation of chondrules is one of the funda¬ 
mental questions of meteoritics and planet formation. 

Chondrules are the 0.1---1 mm-size once-molten sili¬ 
cate spherules found abundantly in primitive meteorites 
known as chondrites (Jones, Grossman & Rubin 2005; 

Sears 2004; Davis et al. 2014). Most of these chondrites 
in fact consist predominantly of these chondrules, so the 
melt-producing events that created them must have been 
extremely common during the first few million years of 
the solar system. Yet there is confusing and conflicting 
evidence as to what these events might have been. Boss 
(1996) published an overview of the status of the discus¬ 
sion at that time, though significant new developments 
have occurred since then. 

Many theories have been put forward over the last half 
a century. One theory involves impact melt sprays. Put 
forward by Urey (1953) and refined by Kieffer (1975) this 
model states that high-velocity impacts (> 3- • • 5 km/s) 
could lead to sprays of impact melt that produce droplets 
which solidify into chondrules. While these required colli¬ 
sion velocities are high, a small fraction of impacts might 
acquire such velocities (Bottke et al. 1994). Perhaps the 
so-called “Grand Tack”-scenario of Walsh et al. (2011), 
in which Jupiter temporarily entered the asteroid belt re¬ 
gion before migrating back outward, might produce such 
high velocities. However, Taylor et al. (1983) put forward 
a number of arguments against the impact-melt scenario, 
some of which were based the differences between chon¬ 
drites and lunar impact regolith. 

Zook (1980) suggested that if the interiors of the collid¬ 
ing bodies are already in a molten state due to 26 Al de¬ 
cay, then the required impact velocities to create sprays 


of melt are much lower, and thus more consistent with ex¬ 
pectations of the average relative velocities between plan¬ 
etesimals. Sanders & Taylor (2005), Hevey & Sanders 
(2006) and Sanders & Scott (2012) follow up on this idea 
of pre-molten impactors and back it up with models and 
meteoritic evidence. Sanders & Scott (2012) argue that 
this scenario is hard to avoid, given that in the first 2.5 
million years most planetesimals were internally nearly 
fully molten by 26 Al decay heat, and that collisions be¬ 
tween planetesimals were extremely common. Each col¬ 
lision would almost certainly release substantial amounts 
of molten rock into the nebula in the form of sprays of 
lava droplets, and it is natural to assume that these may 
be chondrules. They argue that the near-solar composi¬ 
tion of chondrules can be explained by the vigorous con¬ 
vection inside the molten planetesimals that may slow 
down iron/nickel-core formation, thus keeping the melt 
solar. Some degree of differentiation would then in fact 
explain the low iron-content L and LL chondrites. 

Recently, Asphaug et al. (2011) performed Smooth 
Particle Hydrodynamics simulations for such impacts to 
demonstrate the dynamics of this scenario. In addition 
they proposed a simple way to calculate the melt droplet 
size by equating the released enthalpy after the collision 
to the surface energy of the droplets. 

There are many alternative scenarios proposed in the 
literature. Perhaps the most popular model is the flash 
heating by nebular shocks (Hood & Horanyi, 1991). De- 
sch & Connolly (2002) and Ciesla & Hood (2002) de¬ 
veloped detailed 1-D radiative shock models with dust 
particles and chondrules interacting both radiatively and 
frictionally with the gas. They showed how the radia¬ 
tive shocks exhibit temperature spikes of mere tens of 


2 


Dullemond, Stammler, Johansen 


seconds (with cooling rates > 10 4 K/hr) that would be 
good candidates for chondrule-forming events. After the 
main shock temperature spike their model exhibits a fur¬ 
ther cooling at intermediate cooling rate (~50 K/hr), 
which cools the chondrules to sufficiently low tempera¬ 
ture in a sufficiently short amount of time. This model 
seems to produce the flash-heating events required for 
turning “dust bunnies” into chondrules. However, so far 
the shock scenario still has several unresolved issues (e.g., 
Desch et al. 2012; Boley et al. 2013; Stammler & Dulle¬ 
mond 2014). 

Also some issues are raised about whether large scale 
shocks in the optically thick solar nebula can explain the 
time scales involved in chondrule formation (Stammler 
& Dullemond 2014). It is noteworthy that there exist 
several other nebular flash heating models, most notably 
flash heating by nebular lightning (Gibbard et al. 1997) 
and flash heating by energy dissipation in current sheets 
forming in MHD turbulence (Hubbard et al. 2012). 

One important constraint on chondrule formation 
models is that these heating events were very short-lived 
compared to any other nebular time scales. By compar¬ 
ing textures of chondrule to those obtained in furnace 
experiments it can be inferred that chondrules must have 
cooled from the liquidus temperature down to below the 
solidus temperature in a matter of hours (Hewins 1983; 
Hewins & Connolly 1996; see also references in Morris & 
Desch 2010 and the excellent review paper of Desch et 
al. 2012). In other words: the chondrule forming events 
must have been flash-events. On the other hand, un¬ 
der optically thin conditions a molten chondrule would 
radiatively cool in about a second, which would be too 
fast. 

Another constraint that a chondrule formation model 
must fulfill is the retention of volatile elements such as 
Fe, Mg, Si, Na and K. Chondrules are not observed 
to have low abundances of these elements. In particu¬ 
lar for the highly volatile elements Na and K this is a 
puzzle, because any heating event that heats a (proto- 
) chondrule above ~ 1700 K and keeps these chondrules 
above that temperature for more than a few minutes 
will cause most of the Na and K to evaporate out of 
the chondrule. Alexander et al. (2008) argue that this 
means that chondrules must have formed in regions that 
are extremely dense in solids, much more so that typi¬ 
cal dust concentration mechanisms in the protoplanetary 
disk can achieve. According to their calculations such 
regions are so dense in solids that they must be self- 
gravitating. Morris et al. (2012) instead propose that 
chondrules formed behind the bow shock of a fast moving 
planetary embryo, and that the embryo’s atmosphere is 
rich in such volatile elements caused by outgassing from 
the embryo’s interior, thus providing the necessary vapor 
pressure to keep the chondrules volatile-rich. 

In this paper we will focus on the impact splash hy¬ 
pothesis, either the low-velocity pre-molten planetesimal 
version or the high-velocity inrpact-melt version. We as¬ 
sume that after the collision between two planetesimals 
a cloud of molten droplets of magma was released. Since 
soon after the ejection of this cloud of magma droplets 
the internal pressure of the cloud would have dropped 
to near-zero (there is, to good approximation, only vac¬ 
uum between the droplets), the cloud would simply ex¬ 
pand ballistically. Initially the density of the cloud is 


so high that the cloud is completely optically thick and 
no radiation can escape from its interior. There is also 
no adiabatic cooling because of the lack of pressure. So 
the temperature of the magma droplets will initially stay 
roughly constant in time. As the cloud expands, how¬ 
ever, the optical depth drops and eventually the cloud 
will start to cool radiatively. During the early stages the 
density is very high and the evaporation of volatile ele¬ 
ments will be saturated (as was argued by Alexander et 
al. 2008). Once the temperature drops below the solidus, 
volatiles can no longer escape. 

In this first paper we intend to compute how the tem¬ 
perature behaves as a function of time after the impact. 
To do this we set up a simple model: that of a ho¬ 
mologously ballistically expanding homogeneous spher¬ 
ical cloud of lava droplets. We compute the temperature 
of the chondrules as a function of time t and comoving ra¬ 
dial location inside the cloud r/r c i ou d by solving the time- 
dependent radiative transfer equation. We also present 
an analytic approximate solution. While the spherical 
cloud model is not an accurate model of the complex 
shape of an impact splash, it can be regarded as a model 
of part of the impact splash. As such we believe that 
it will give reasonable estimates of the radiative cool¬ 
ing behavior of the impact splash as a function of model 
parameters such as the total mass of the cloud and its ex¬ 
pansion velocity. The total mass of the cloud tells some¬ 
thing about the masses of the two colliding bodies while 
the expansion velocity of the cloud tells us something 
about the impact velocity. 

2. EXPANDING CLOUD MODEL 

When lava droplets are produced in a collision between 
two planetesimals, they will disperse away from the im¬ 
pact site in a ballistic way. Let us call this the “impact 
splash”. If the impact velocity is larger than the escape 
velocity of the two colliding planetesimals this ballistic 
(pressureless) expansion will be linear (i.e. the velocity 
of expansion will not change with time). We regard the 
impact splash as an expanding cloud of lava droplets 
(chondrules-to-be). This cloud is not necessarily cen¬ 
tered on the impact site; it can also move away from the 
impact site. The radius of the cloud increases linearly 
with time t and thus the density of the cloud will drop 
as 1/t 3 . Depending on the complexity of the geometry of 
the impact splash (see e.g. Asphaug et al. 2011) we can 
also consider the splash to be multiple smaller clouds. 
In either case, each of these clouds will expand linearly. 
The scenario is pictographically shown in Fig. [l] 

If the impact is at low velocity, then the gravitational 
pull of the surviving (or merged) body can cause some 
(or most) of the material to re-accrete. In this case the 
1 /t 3 expansion will cease and turn into recompression. 
Initially, however, the linear expansion and 1/t 3 density 
drop will still be a good description. 

To compute the temperature of the droplets in these 
clouds as a function of time after impact we must solve 
the time-dependent equation of radiative transfer. This 
is a very non-local problem, because a cooling droplet 
at location x\ radiates away some heat into the form of 
infrared radiation, which can then be absorbed by an¬ 
other droplet at location x 2 7 ^ x\ which is then heated. 
The droplets are thus radiatively coupled over large dis¬ 
tances. Solving this problem in 3-D for complex cloud 
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Pre-collision After collision: Some parts survive. 

Debris clouds (e.g. clouds of molten 
lava droplets) get ejected. 


The debris clouds (here modeled as 
spherical clouds) expand linearly as 
they move away from the impact site. 


Fig. 1.— Pictographic representation of the model. Left: Two planetesimals approach each other and are about to collide. Middle: After 
the collision some parts of the planetesimals may survive, but some parts are destroyed and dispersed in a cloud of debris. We model this 
dispersing cloud as a set of spherical clouds. If the debris consists of melt, the cloud will consist of molten droplets (chondrules). Right: 
As the clouds move away from the point of impact they expand linearly. The average distance between the chondrules increases linearly 
with time. 


geometries is somewhat challenging. 

We believe, however, that much can already be learned 
from a simple spherical cloud model, for which the radia¬ 
tive transfer problem can be solved to high precision and 
reliability with the 1-D tangent-ray variable eddington 
factor method. We will describe the method in Section 

SI 

The three main parameters of the model are: 


the cloud is constant and outside of the cloud is zero. 
The density of the cloud as a function of time is then 


( 3 A/,,] mid 

/Odoud(r,t) = 47rfl Lud( 4 ) 


for r < i?cioud(t) / p. \ 

for r > i?doud(t) 


From here on we write only p c ioud(t) instead of 
Pcioud(^,t)- The time-dependence of p c ioud(t) is 


T/ c lo U d , To , ^exp (1) 

where M c i ou d is the total mass worth of chondrules (lava 
droplets) in our spherical cloud, To is the initial tempera¬ 
ture of the chondrules and v exp is the expansion velocity 
of the outer radius of the cloud. To get a feeling for the 
numbers it is more convenient to express the mass of the 
cloud in terms of the radius i? me it,o of a ball of magma 
of the same mass M c i ou( j: 

ALcloud = -^-fchonT melt o ( 2 ) 

where £ c hon is the material density of the chondrule 
droplets and thus, by definition, the material density of 
the hypothetical ball of magma. For a linearly expanding 
cloud the radius of the cloud is 

Tlcloud(t) = V exp t (3) 

where the time t is the time since the impact and the 
ejection of the cloud of melt droplets, assuming perfectly 
ballistic (pressureless) expansion and of course assum¬ 
ing large enough t that i? c ioud(t) i? m eit,o so that the 
droplets are clearly separated from each other. The ve¬ 
locity profile inside the cloud is 

v(r,t)=v ex p r (4) 

-^cloud ) 

which is of course only valid for r < R c \oud(t)- 
Let us define a coordinate r centered on the center of 
the spherical cloud. We assume that the density within 


Pcloud(t) 


3M cloud 1 
4tt?4p t 3 


( 6 ) 


Let us consider chondrules (lava droplets) of radius a c hon 
and material density £ c hon- We take £ c hon = 3.3 g/cm 3 
for our model. The mass of a chondrule is then 


_ 4-71- 3 

TTlchon — “^“Schon&chon 


(7) 


The number density of chondrules is then 


n c hon(t) = (8) 

^chon 


The geometric cross section of a chondrule is i ra(: hon - Let 
us assume that the chondrule has zero albedo. Then 
the absorption cross section equals the geometric cross 
section. The opacity (i.e. the cross section per gram) is 
then 


7ra chon 

47r£ c honachon/3 


3 1 


4 T:[ion Cchon 


(9) 


We will assume that the opacity of any possible vapor 
between the chondrules will be low compared to that of 
the chondrules so that it can be ignored. The optical 
depth from the center of the cloud to the edge is 

3 Mk 1 

T(t ) = Pcloud (t ) K-Rcloud (t) = 

47r ^exp ^ 


The rate of cooling will strongly depend on this optical 
depth. 
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Model 


^exp 

To 

^cool 

'T’cool 

FI 

1 km 

100 m/s 

2000 K 

27 min 

92 

F2 

0.1 km 

1000 m/s 

2000 K 

16 sec 

9.2 

F3 

10 km 

100 m/s 

2000 K 

7 h 

366 

F4 

0.01 km 

1000 m/s 

2000 K 

1 s 

2.3 


TABLE 1 

The model parameters (first three columns) for the main 
(“fiducial”) models presented in this paper. The cloud 

MASS Mcloud CAN BE DERIVED FROM M cloud = ^ £ c hon , 0 

(cf. Eq.EJ. The last two columns are computed from 

THESE PARAMETERS: THE COOLING TIME t coo i (THE TIME AFTER 
THE IMPACT WHEN THE CLOUD STARTS TO COOL, SEE EQ. ITT)! AND 
THE OPTICAL DEPTH AT THAT TIME T coo i = r{t = t cool ) (SEE 

EO,l20ll. 


In this paper we will present our results based on a 
set of fiducial models as well as parameter scans. The 
parameters of the fiducial models are listed in Table |T] 


3. ANALYTIC ESTIMATE OF THE COOLING BEHAVIOR 
OF THE CHONDRULE CLOUD 

As the cloud expands it can start to radiate away en¬ 
ergy. A fully fledged time-dependent radiative cooling 
computation will yield temperature as a function of time 
and space: T(r,t). We will compute this in Section 0] 

A first estimate of the cooling behavior of the cloud 
can already be made with pen and paper. The cloud will 
initially be optically thick: r 1. The luminosity of the 
cloud is then 

L = AnR 2 cloud aT 4 s (11) 

where a is the Stefan-Boltzmann constant and T e g is 
the effective surface temperature at r = R c i ou d- In the 
above equation, and from here onward, we omit the (t) 
for aesthetic reasons, but the time-dependence is still as¬ 
sumed. During the cooling phase the effective surface 
temperature will be lower than the central temperature. 
A commonly used estimation of this effect is: 

L ~ 47ri? c 2 loud aT 4 l (12) 

T 

where T is now the central temperature and 1/r is the 
correction factor to account for the lower temperature 
of the surface T e g, and is based on radiative diffusion 
theory which states that T ~ T X l 4 T e & (for r 1). To 
compute the change of the temperature as a result of the 
radiative loss given by Eq. m we must use the heat 
capacity formula. The total thermal energy stored in the 
chondrule cloud (approximating it, for the moment, as 
an isothermal cloud) is: 

E = M cloud c m T (13) 


where c m is the mass-weighted specific heat of the lava 
droplet. A value of c m = 10 7 erg g _1 K _1 is a reasonable 
value which we will adopt here. The radiative loss time 
scale can now be defined as 


^rad(^) — 


E 

T 


A/cloud Cml 


47Ti? 


2 

cloud 


crT 3 


(14) 


becomes smaller. It is therefore to be expected that at 
early times the temperature of the chondrules remains 
constant. Note that since the cloud consists of liquid 
drops that stay at a constant volume and move away 
from each other, there is no adiabatic cooling involved 
here. Any potential vapor may adiabatically cool, but 
we will ignore this effect in this paper, assuming that the 
total mass in vapor is always small compared to the mass 
in droplets. 

At some point in time, however, t rad becomes smaller 
than t and the chondrules start to cool. Based on the 
results of the true radiative transfer calculations of Sec¬ 
tion 0] it turns out that near the center of the cloud the 
transition from the initial constant temperature phase to 
the temperature-decline phase occurs roughly at a time 
tcooi defined by 


1 rad (tc 


= 5 t c 


>ol 


Inserting Eas. (U4ll2ll3llOD into Eq. (fl5l) yields 


/I 3 t\I^ loud c m h- \ 
V5(4tt) 2 u 4 xp <rT 3 ) 


(15) 


(16) 


Note that in Eq. (fT5l) the factor of 5 is purely empirical. 
The estimate we make here is essentially a dimensional 
analysis in which proportionality factors have to be cal¬ 
ibrated against more exact calculations (in our case the 
full radiative transfer calculation). 

Roughly for t < t coo l the temperature stays constant 
while for t > t CO oi the temperature drops with time. 
From now on we shall define the t coo i as the one cal¬ 
culated for the center of the cloud. Near the surface the 
cooling sets in earlier. 

The temperature decline with time can also be es¬ 
timated, at least up until the point where the cloud 
becomes optically thin, after which our approximation 
breaks down. The way to make this estimate is to solve 
the cooling equation 


dT(t) _ T(t) 

dt tcooi (t) 

where t coo i(t) is the cooling time given by Eq. (fTfil) but 
with To replaced by T(t). The reasoning behind this is 
that if we would instead use t lad (t) (which might look 
more reasonable at first sight) we will quickly cool to 
temperatures for which condition f rad < t is again no 
longer fulfilled. This is because t rad oc T~ 3 and thus very 
rapidly rises with declining temperature. If we insert 
Eq. (fTTH) with T 0 replaced by T into Eq. ijTTll we obtain 
a differential equation that can be solved by separation 
of variables. The solution is 


T(t > tcooi) = T 0 


3 t 


5 tcooi 



(18) 


where here f C ooi is again the original one with T d . We 
assume that for t < t c 00 i we have: 


T(t < tcooi ) = T 0 (19) 


This time scale varies with time: it is very long at early 
times because the cloud is then still extremely optically 
thick. At late times the cloud is optically thinner and 
the radiation can more freely escape (Eq. fl2l) and f ra d 


The solution Eqs. (118I19D holds true for the central tem¬ 
perature. But most of the mass of a homogeneous 
sphere resides in the outer parts. For the temperature 
at r = 0.81?cioud and r = 0.9i? c i ou d (roughly the radii 
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time/t coo . 


Fig. 2.— The temperature evolution of the ex pandin g cloud 
model according to the analytical estimate of Eqs. 1181191 . shown 
in black. In grey are, for comparison, the full radiative transfer so¬ 
lutions from Section [|] for the fiducial model FI (see Table [TJ). The 
temperatures are shown at three positions in the cloud: the center, 
at 80% of the radius and at 90% of the radius ( near the surface of 
the cloud). The time scale is scaled to t coo l fEq. H6D . Th e analytic 
solutions are only valid as long as r coo i 1 (cf. Eq. I20D . 

of half mass and of 75% mass respectively) the best fit¬ 
ting solution to the real solutions of Section 0] are similar 
to Eq. m but with the term 2/5 replaced by 3/5 and 
3.8/5 respectively. These analytic estimates of the tem¬ 
perature as a function of time are plotted in Fig. [2j As 
one can see, although the time after the impact at which 
the cooling stars is different depending on whether you 
look at the center or at the edge of the cloud, the typical 
cooling rate after the cooling begins is similar throughout 
the cloud. 

It is important to remind ourselves that the cooling 
normally starts well before the cloud becomes optically 
thin. Optically thin cooling would be extremely fast: a 
single chondrule would cool within about 1 second. The 
protracted cooling over a time scale of hours occurs be¬ 
cause the high optical depth (even after t = t c 00 i) acts 
as a kind of blanket that keeps the chondrules warm, al¬ 
beit a blanket that is thinning over time as the optical 
depth drops. To get a feeling for this we can calculate 
the optical depth at the time t = t coo i from Eq. m- 


o3/5 5 2/5 M 1/5 3/5 2/5 T 6/5 

^ \ — ° ° iM cloud^ a / on \ 

Tcooi - T (tcooi) - 2/5 2/5 

(47T y/°V e kpCrn 


So Eqs. (11SI19I) . which are based on the assumption of 
the “blanket effect” of high optical depth, are expected 
to be valid for the case t coo i 1. In table Q] one can see 
that T CO oi is well above 1 for all models, except the most 
extreme model F4 which has anyway a much too small 
cooling time to be consistent with chondrule textures. 

The analytic model of Eqs. (1181191) tells us that for t > 
t C oo\ the temperature drops off suddenly on a typical time 
scale of tcooi- The temperature decay rate at t = t coo i is: 


T(t-fcooi)- dt 


To 


t -^cool 


Dol 


( 21 ) 


This is plotted, together with t coo i, in Fig. [3] It is also 
useful to define another kind of cooling time scale ti 4 oo 


which is the time between the start of the cooling at 
T = To and the time when the temperature has dropped 
below the solidus temperature 1400 K. To first order we 
can write 


T 0 - 1400/v 

tl400 — —:- 

\T(t=tcool)\ 



1400AT \ 

To J 


tcooi 


( 22 ) 


The quantities T and tuoo can be compared to the con¬ 
straints coming from the analysis of textures of chon¬ 
drules, which put time limits on the cooling process of 
order of hours, or more precisely: cooling rates in the 
range 10 • ■ ■ ~ 3000 K/hr (see Morris & Desch 2010 and 
Desch et al. 2012 and references therein, though see also 
Miura & Yamamoto 2014 for a different view based on 
theoretical modeling of crystallization), indicated by the 
grey area in Fig. [3j According to Eq. (fl6l) i CO oi goes as 

ilffkmd and as v e J^ 5 - There is not very much freedom 
of choice of To: it must lie somewhere between 1770 and 
2120 K (see Morris & Desch 2010 and references therein). 
The n does not have too much wiggle room either: the 
radii of chondrules are known. This means that con¬ 
straints on f coo i from textural analysis directly set lim¬ 
its on the ratio M c i ouc j/v 2 xp or equivalently on the ratio 

-^meit o/^exp- for example, we choose a c hon = 0.03 
cm and Tq = 2000 K we obtain n = 7.6 cm 2 /g. If 
we require, for instance, the cooling time scale to be 
Aooi = 10 minutes, then we find that M c i ou d/u 2 xp ~ 10' 
g s 2 cm -2 . This would be fulfilled e.g. for r) exp = 1 km/s 
and Mcioud = 10 17 g (a mass corresponding to a ball of 
magma of roughly i? me it,o =2 km radius). It would also 
be fulfilled by a smaller mass and smaller velocity: e.g. 
for u exp = 10 m/s and M cloud = 10 13 g (i? me it,o =0-1 km 
radius). Such small mass could either mean that upon 
impact only a small fraction of the debris is in the form 
of chondrule droplets, or it could simply mean that the 
largest closed unit of the droplet splash (subcloud) is so 
small, but many such expanding cloudlets of chondrules 
are ejected. The elongated streams of debris found in 
Asphaug et al. (2011) could be regarded as a string of 
such smaller mass cloudlets. 

This is shown more directly in Fig. [4] which shows the 
parameter space of the model, with the grey area again 
showing the models which give cooling rates between 10 
and 3000 K/hr. 


4. TIME-DEPENDENT RADIATIVE TRANSFER 

So far we have only made an estimate of the cooling 
behavior of the cloud of liquid chondrules. Let us now 
calculate the temperature profile with a full treatment of 
time-dependent radiative transfer. The time-dependence 
comes in only due to the time it takes the lava droplets to 
convert their heat into radiation, not in the form of light- 
travel time (which typically only plays a role for media 
at temperatures above 10 5 K). This can be easily verified 
by comparing the radiative energy density aT 4 to the 
material energy density pc m T , which in our case always 
satisfies aT 4 -C pc m T. Therefore the radiative transfer 
equation itself is stationary, and the time-dependence is 
in the equation for heating/cooling of the chondrules. 


4.1. Equations of time-dependent radiative transfer 
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Fig. 3.— The cooling time t CO ol (Eq. m as a function of param¬ 
eters of the model for To = 2000 K and a c hon = 0.03 cm. Note 
that the time t o co ol down from 2000 K to 1400 K (£ 1400 ) is, ac¬ 
cording to Eq. l22l I . about four times smaller than t coo l- On the 
right axis the corresponding cooling rate \dT/dt\ in K/hr is plotted. 
The grey area is the observed range of chondrule cooling rates (see 
main text). Note that the values of £ C ool are (by definition) for the 
center of the cloud. The regions of the cloud more closely to the 
surface will start to cool earlier, but the cooling rate is almost the 
same. 



Fig. 4.— The parameter space of the model. The grey area 
shows where the solutions have cooling rates in the range observed 
for chondrules. In the bottom-right the region is shown where 
T cool < 10, where the analytic solution becomes invalid. 

The formal radiative transfer equation is (see e.g. Mi- 
halas & Mihalas 1999): 

n ■ VI(x,n) = pn(B(x) — I(x,n)) (23) 

where I is the frequency-integrated mean intensity in erg 
s^ 1 cm -2 ster -1 , x is the position vector in space, n is 
the unit vector of direction of the radiation, and finally 
B the frequency-integrated Planck function given by 

B(T) = -T 4 (24) 

7T 

We will use frequency-integrated quantities because the 
opacity of a chondrule is expected to be roughly equal 


to the geometric cross section, independent of frequency, 
because the chondrules are much larger than the wave¬ 
length of infrared radiation. In this case the frequency- 
dependence of the radiative transfer does not need to 
be considered. Eq. (l23l) is called the “formal transfer 
equation”, and describes the transport of radiation along 
rays of direction n. It is easy to integrate if the value of 
B{x) = B(T(x)) is known everywhere, but since we want 
to calculate the temperature T(x), this function is part 
of the thing we want to solve. 

In our 1-D spherically symmetric setting Eq. (l23l) can 
be written as 


: . dj /4 r ) 

dr 


1 ~ jJ 2 dip (r ) 
r dp 


a(B(T(r))-/ M (r)) (25) 


where we shortened 


pK =: a (26) 

Here p = cos (0) where 9 is the angle between the ray 
along which the radiation is followed and the radially out¬ 
ward pointing unit vector. The Eddington factor method 
of solving this equation relies on the definition of the first 
three angular moments of the intensity at each location 
r: 


J i r ) = \ J I ni r ) d » (27) 

H ( r ) = \ J (28) 

K i r ) = \ J i I^{r)p 2 dp (29) 


Here J is called the mean intensity and H is the flux 
divided by 47 t. By integrating Eq. (1251) over ^dp after 
multiplying it by 1 and by p respectively one obtains the 
first two moment equations: 


1 d{r 2 H) 


dK 

dr 


2 dr 

3/-1 


= a(B(T) — J) 


J = —aH 


(30) 

(31) 


These are two equations with three unknowns. To close 
this set of equations we write 


I<(r) = f(r)J{r) (32) 

where / is called the “Eddington factor”. If /(r) is 
known, then the two moment equations are in closed 
form and can be solved for J(r) and H(r). From the 
moment equations alone we cannot know what f(r) is 
for all r, but for certain limiting cases we know their 
values: For optically thick media / = 1/3. Outside of 
the cloud for r —>• oo we will get / —> 1, which is the 
free-streaming limit. However, for general r and for gen¬ 
eral values of the optical depth we must employ another 
method of calculating f(r). We do this by integrating, 
for the current value of T(r) (and thus B(T(r))), the 
formal transfer equation Eq. (EH) using the “tangent ray 
method” (see e.g. Mihalas & Mihalas 1999). By inte¬ 
grating the resulting intensities over angle we calculate 
the current estimates of J(r) and K(r). Let us call these 
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Jfte(r) and JtTfte(r) (where “fte” stands for “formal trans¬ 
fer equation”). Our current estimate of /(r) is then 


f ( r ) 


Khe(r) 


(33) 


obtain J" . The problem with this scheme is that the 
required time steps can become very small due to nu¬ 
merical stiffness. A more robust method is to integrate 
the complete system Eos. (1351361) using an implicit inte¬ 
gration scheme. This is discussed in appendix [A] 


This is then the /(r) function we stick into Eq. (1521) . so 
that the moment equations can be solved for J(r) and 
H (r). We also need to impose boundary conditions at the 
inner and outer edge. At the inner edge we take H = 0 
(zero flux). At the outer edge we could set H = J, which 
is valid if the outer edge of our computational domain is 
at r —>■ oo. For finite outer radius we set instead 


H(r 0 ut ) 


77 ft e ('/'out) 

d fte (r'out) 


j (I’out) 


(34) 


Note that we can choose to set the outer edge of our 
computational domain r ou t at r out = i? c i ou d, but we can 
also set it at ?' ou t > -Rdoud, as long as we properly set the 
boundary condition at r = r ou t according to Eq. (El. 

We can now combine the two moment equations into 
the following form: 

lL±{ r l WJ) , 3 f-l j 

a r 2 dr \ a dr r 


J = J ^T 4 (35) 


This can be solved, together with the boundary condi¬ 
tions, using a matrix equation. The details of this are 
discussed in Appendix [Al 

Next we must compute how the temperature T(r) re¬ 
acts to the radiation field, or in other words, how T(r) 
radiatively cools with time. The energy equation is 

dT ( a 4 \ 

p Cm —=4^(j--T) (36) 

where p is given by Eq. ©• The d/dt operator is the 
comoving derivative, i.e. the time derivative is computed 
for a given chondrule. Eq. (l36l) is a local equation at 
each location. The non-locality of the heating/cooling 
is only through the non-locality of the equation for J(r) 
(Eq.SSJ. 

The way we solve these equations numerically is by 
setting up a radial grid of N gridpoints ri (with 1 < 
i < N) upon which we define the mean intensity Ji and 
the temperature X). Since the cloud is expanding we let 
the gridpoints move along with the material (Lagrangian 
approach) so that the location of grid point i at time t is 

r i(t) = ViVexpt (37) 

where r/ is the dimensionless radial coordinate 

T T 

r, = - -— =- (38) 

A c loud(t) Uexpt 

The values of ip do not change with time. Each chondrule 
stays at constant rj, and so the comoving time deriva¬ 
tive d/dt used in e.g. Eq. (13^1) becomes simply the time 
derivative of that quantity at a given r/i grid point. 

In principle one could imagine a simple time-dependent 
numerical integration method for the combination of 
Eos. ((351361) . We know J" and T” at time t = t n , and 
we can solve for Tf/ +l by taking an explicit Euler inte¬ 
gration step of Eq. (IM1) . Then we can solve the moment 
equations (because now we have B/ + = B(T/f + )) and 


4.2. Results of the time-dependent radiative transfer 

models 

The results of the time-dependent radiative transfer 
calculations are shown in Fig. [5] These results confirm 
the estimates made in Section [3] We find that the shape 
of the cooling curves is mostly the same for all models, 
except for the time scaling, which we know from the ana¬ 
lytically expression of t coo i (Eq- ESI) ■ Only for the rather 
extreme (and presumably unrealistic) case such as fidu¬ 
cial model F4 with = 0.01 km and u exp = 1000 

m/s we find from the full radiative transfer model that 
the cooling is slower than the analytic solution. This is 
because in that case the condition that r coo i (Eq. O is 
1 is broken and the cooling time scale will then not be 
determined by the optical depth decline but by the effec¬ 
tiveness by which a chondrule can convert its heat into 
radiation. These cases, however, imply typically time 
scales of seconds rather than hours, and these are any¬ 
way too short to be consistent with chondrule textures. 

It seems that the analytic solutions are fairly accurate 
when t coo i > 10, which is always true for the interesting 
regime of parameter space. 

From Figs. [2] and [5] we see the model systematically 
predicts a period of constant temperature with a dura¬ 
tion similar to the later cooling phase. This is true at 
least near the center of the cloud. Near the surface this 
temperature plateau is shorter, but the cooling time is 
similar. 

In Fig. [6] we show, for model FI, how the temperature 
looks as a function of radial position in the cloud, for 
different times after the collision. Since the cloud is ex¬ 
panding, which would make it difficult to compare the 
temperature profiles at different times, we use the radial 
coordinate relative to the outer cloud radius. One sees 
that the outer regions start cooling earlier because the 
optical depth to the surface is smaller. Also one sees 
that a negative radial temperature gradient is produced 
which causes the radiative diffusive energy transport to 
the surface. 

The cooling curve from our model is significantly dif¬ 
ferent from the ones predicted from shock heating (e.g. 
Morris & Desch 2010, Morris et al. 2012). The shock 
heating models predict a radiative preheating phase, a 
temperature spike with rapid cooling, and depending on 
geometric conditions followed by a slower cooling phase. 
In our model the temperature is high from the start (or 
for high-speed collisions: upon impact), stays relatively 
constant for some time and then drops over a similar 
time scale. 


5. DISCUSSION 

The scenario of splashes of melt droplets originating 
from colliding pre-molten planetesimals has been dis¬ 
cussed at length in several papers by Sanders and cowork¬ 
ers (e.g. Sanders & Scott 2012) and the paper by Asphaug 
et al. (2011). These papers discuss how this scenario 
holds up agains the many meteoritic constraints. We 
will not repeat these arguments here, but refer for those 
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Fig. 5.— The temperature evolution of the fiducial expanding cloud models FI (top-left), F2 (top-right), F3 (bottom-left) and F4 
(bottom-right). See Table [T] for the parameters. The temperatures are shown at three positions in the clou d: t he center, at 80% of the 
radius and at 90% of the radius (near the surface of the cloud). The vertical line shows the time t coo l of Eq. csj. 



r/Rci„ U d(t) 


Fig. 6.— The radial temperature temperature structure of fidu¬ 
cial model FI for different times. The horizontal axis is the radius 
relative to the outer cloud radius. Note that the temperature pro¬ 
file extends beyond the outer radius of the cloud because even if 
that region is empty, a chondrule residing there would have a well- 
defined temperature. For our purposes, however, only the temper¬ 
ature within the cloud is relevant. 

discussions to those papers. Instead, we will focus on the 
discussion of aspects related to our splash cooling model. 


5.1. Splash geometry 

Our model describes a simple spherical expanding 
cloud as a model for (part of) the splash resulting from 
colliding bodies. It is obvious that this is a drastic simpli¬ 
fication. The SPH model of Asphaug et al. (2011) shows 
a much more realistic geometry involving, among other 
complexities, tube-like geometries. A proper 3-D time- 
dependent radiative transfer simulation done in concert 
with the gravito-hydrodynamic simulation of the impact 
splashing is required to know precisely what the cooling 
behavior is. Nevertheless we predict that our result can 
be used to estimate the results of those detailed compu¬ 
tations: the tube-like geometry can be approximated as 
a homologously expanding cylinder, expanding not only 
in the 2 perpendicular directions but also in the longitu¬ 
dinal direction. Its properties are then not too much dif¬ 
ferent from a chain of spherically expanding clouds with 
diameters similar to the diameter of the tube, except 
that these spherical cloud cannot cool in all directions 
(47 t) but only in part of the sky (in the other part it will 
heat its neighbors and be heated by its neighbors). We 
would then expect that the cooling time f coo i will be a 
bit longer but not by a large factor. 

Also, one can expect some parts of the splash to be 
ejected at slow speed compared to the escape speed, 
so that it may fall back to the remaining (unsplaslied) 
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bodies and accrete on them. This evidently breaks the 
assumption of homologous expansion, and as a result 
will also deform any initially spherical sub-cloud into 
a more pancake-like shape. The idea of regarding the 
splash cloud as a collection of spherical sub-clouds is then 
clearly no longer valid, not even in an approximative way. 
In such cases a full 3-D radiation-gravito-hydrodynamic 
simulation of the splash is unavoidable. However, multi¬ 
dimensional radiation-hydrodynamics is an enormously 
challenging problem, in particular for problems involving 
such extremely complex geometries. It is an interesting 
challenge for the future. For getting order-of-magnitude 
estimates we may be forced to stick, for the time be¬ 
ing, with simple models such as the one presented in this 
paper. 

5.2. Chondrule sizes 

According to the model by Asphaug et al. (2011), by 
which the lava drop size is calculated from the initial 
pressure under which the magma was in the pre-molten 
planetesimals before the collision, the planetesimal sizes 
must have been 10 km or larger. Since our cloud of lava 
droplets can be one of a number of such clouds ejected 
from the impact, the i? me it,o of our model cannot be di¬ 
rectly compared to that > 10 km size lower limit from 
Asphaug. But it does give an indication that, if we adopt 
their model of the chondrule sizes, it seems reasonable to 
be looking to values of roughly the order of i? me it,o — 10 
km or larger. To elucidate our reasoning: If we collide 
two fully molten planetesimals of 20 km radius, and we 
model the splash with 10 spherically expanding droplet 
clouds, then we would get -R me it,o — 12 km. In Fig. [3] 
this means we are looking at the region of the diagram 
above the dotted line. 

There is a caveat, however: What tells us that it is just 
10 clouds, not a 1000? This brings us to the complicated 
issue of the actual geometry of the splash. If we “de¬ 
compose” the splash geometry into many small spherical 
clouds we must assure that these clouds can radiatively 
cool to the outside. In other words: we cannot consider 
a big spherical cloud as consisting of many smaller sub¬ 
clouds, because they will irradiate each other and not 
cool much. But if the cloud geometry is a long extended 
streak, it can be considered as consisting of a chain of 
spherical subclouds. We are here confronted with the 
limitation of our spherical expanding cloud model. 

5.3. Partial pre-melting 

The reason why pre-melting seems to be necessary for 
the impact splash scenario to work is that very high speed 
collisions (> 5 km/s) needed to produce melt from cold 
planetesimals are presumably rare. However, pre-melting 
may lead to differentiation of the planetesimals, causing 
the iron and nickel to sink into the core and produc¬ 
ing a basaltic mantle. Chondrules formed from such ob¬ 
jects would be very non-solar in composition. Although 
Sanders & Scott (2012) argue that convection reduces 
the efficiency of the differentiation and that some amount 
of differentiation is in fact consistent with e.g. L and LL 
chondrites, it seems that this issue is not yet conclusively 
solved and more detailed and quantitative modeling of 
the differentiation process in these molten planetesimals 
may need to be done. 


However, a middle way is if planetesimals were pre¬ 
heated to elevated temperatures below the melting point. 
This is to be expected in particular for later generation 
planetesimals (i.e. planetesimals formed at a time when 
a substantial fraction of the 26 Al has already decayed). 
Lower velocity collisions might then be sufficient to add 
the remaining energy needed to produce melt. The prob¬ 
lem with this scenario is that it requires fine-tuning: The 
planetesimal must have formed at a time when the re¬ 
maining abundance of 26 Al is low enough not to melt 
the planetesimal but high enough to heat it to tempera¬ 
tures only slightly below the melting temperature. This 
might be too much coincidence. 

On the other hand, we know that some planetesimals 
were strongly differentiated: the parent bodies of iron 
and basaltic meteorites. If chondrules formed from col¬ 
lisions between planetesimals, one would expect some 
chondrules to have been formed from collisions with such 
differentiated planetesimals. Some chondrules should 
have compositions that are very iron poor, presumably 
more iron-poor than LL-chondrites. This raises the ques¬ 
tion where are these? One possible answer is that chon¬ 
drules produced in these very early phases (< 1 million 
years after CAI formation) were either accreted into the 
sun or reincorporated into other bodies, where they were 
again molten. But if this process of elimination is so ef¬ 
ficient, why did some CAIs survive and get incorporated 
into the (later) chodrites? This shows that while the im¬ 
pact splashing model is appealing, there is more work to 
be done to answer such questions. 

5.4. Compound chondrules 

In our homologously expanding cloud model the melt 
droplets move away from their nearest neighbor at ex¬ 
tremely low velocity (millimeters per minute). Small in¬ 
homogeneities (which are of course unavoidable in such a 
messy process as an impact) would lead to some droplets 
actually colliding. When they do so during the initial 
constant temperature phase, they presumably coalesce 
and form a larger melt droplet. However, when they col¬ 
lide during the cooling phase, they may have become cool 
enough to be plastic but not liquid anymore. They will 
then produce a compound chondrule. 

It would be useful if we could quantify the percentage 
of chondrules that become compound in this way. How¬ 
ever, at present we see no way how we can estimate the 
random relative velocities between adjacent chondrules 
in the expanding clouds. 

Many compound chondrules, however, clearly show one 
chondrule having remained rigid, while the other having 
been indented. So clearly they were of different temper¬ 
ature. This is difficult to achieve in our simple spherical 
expanding cloud model. For this it will be necessary to 
have other debris (e.g. another cloud-fragment emerging 
from the impact site) to interdisperse with our cloud of 
droplet, so as to have chondrules from different thermal 
environments to mix. And this must happen on a time 
scale shorter than the time scale of thermal equilibration 
between the chondrules from these different origins. To 
make this more quantitative would be, however, a chal¬ 
lenge. It would likely require us to contemplate complex 
splash geometries. 

Metzler (2012) found that some chondrites in fact con¬ 
tain entire large scale clusters of compound chondrules. 
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A possible explanation could be that some of the im¬ 
pact debris may have fallen back onto the remainder of 
the biggest of the two colliding planetesimals. If this 
happens quickly enough, so that this re-accretion occurs 
before the chondrules have cooled below the temperature 
at which the chondrules are still sticky, then it is reason¬ 
able to assume that such clusters may form as chondrules 
simply fall onto each other. 

5.5. What about relict grains and rims? 

Another problem with our simple spherically symmet¬ 
ric expanding plume model is that it does not account for 
any possible “pollution” of the melt droplets with other 
(non-molten) debris. Some chondrules contain relict 
grains (Jones 1996), which must have been inserted while 
the chondrule was still liquid, but the chondrule must 
have cooled very soon thereafer. As in the case with the 
compound chondrules (Section 15.41) this would require 
debris from a different part of the impact debris cloud to 
interdisperse with our cloud of droplets. This foreign de¬ 
bris must have originated from non-molten parts of the 
original planetesimals, e.g. the regolith or pristine dust 
covering the crust of the original planetesimals. Only 
those grains from that non-molten debris that entered 
late enough can enter the chondrules shortly before they 
cool below the solidus, and thus survive. Like with the 
compound chondrules, however, getting quantitative es¬ 
timates of this process is challenging. 

5.6. Producing chondritic parent bodies 

The impact model may provide a plausible scenario for 
the formation of chondrules, but how do these chondrules 
(and matrix) accrete to form another planetesimal (the 
parent body of the later chondritic meteorite)? A nice 
aspect of the low-velocity collision of pre-molten plan¬ 
etesimals model is that much of the debris may, in fact, 
reaccrete onto one of the surviving original planetesimals 
(or what remains of it). The idea is that due to the 
high gas density in the protoplanetary disk the stochas¬ 
tic planetesimal motions are of low velocity so that when 
they collide, they collide essentially at velocities that are 
only just above their mutual escape speed. Since much 
of the energy is dissipated in the inelastic impact, most 
of the debris will not have enough kinetic energy to es¬ 
cape the system and will (after perhaps a flyby or two) 
eventually accrete onto a single body or a binary body. 
However, there is evidence that inside a single chondrite 
there exist chondrules of different ages (Villeneuve et 
al. 2009). And some chondrules experienced multiple 
heating events. Somehow chondrules produced in multi¬ 
ple chondrule-forming impacts must have mixed together 
and form a single chunk of rock that later is observed as 
a chondritic meteorite. The immediate re-accretion sce¬ 
nario may make this somewhat difficult because once the 
collision releases most of the magma, and much of this is 
reaccreted as chondrules (for the molten part) and ma¬ 
trix (for the non-molten part) the remaining body will 
no longer be a sphere of magma. 

If, however, the collisions would be high-speed enough 
that much of the debris escapes, the newly formed chon¬ 
drules will disperse into the protoplanetary disk. There 
they mix with the pool of other chondrules of various 
ages, and then later accrete into a chondrite parent body. 


An important question in this scenario is how aerody¬ 
namic drag and the resulting radial drift affects this: 
will the chondrules stay long enough in the asteroid-belt- 
region before they dift away or not? Jacquet et al. (2012) 
study this question in detail. 

Perhaps more likely is a combination of the two. In 
particular the cluster chondrites of Metzler (2012, see 
Section 15.41) can be easily understood in terms of the 
reaccretion scenario while the multiple ages seems more 
to point toward a dispersion and later accretion scenario. 

5.7. Volatile elements 

One of the main characteristics of most chondrules is 
their “normal” abundance of volatile elements such as 
Na, K etc. If a chondrule would stay at a high tempera¬ 
ture (e.g. 2000K) for more than a few minutes, these el¬ 
ements would have evaporated out of the chondrule and 
leave a volatile-depleted chondrule behind. We believe 
that the impact splash scenario may naturally resolve 
this issue, because during the hot phases of the expand¬ 
ing impact splash the density of the cloud is so high that 
the vapor quickly reaches the saturation pressure and 
evaporation stops. However, this idea must be verified 
by actual modeling: Will the evaporation indeed satu¬ 
rate? Will the temperature drop below the evaporation 
temperature before the density of the expanding cloud 
becomes too low to saturate the pressure? Answering 
these questions with explicit calculations is the topic of 
paper II in this series. 

5.8. What kind of chondrules qualify for this scenario? 

So far we have not made any distinction between differ¬ 
ent kinds of chondrules or chondrites. Given that chon¬ 
drule properties vary considerably between chondritic 
classes, the natural question arises: do all chondrules 
form through the same mechanism or are different chon¬ 
drule formation scenarios responsible for different kinds 
of chondrules? And for which kind, if at all, could the 
impact splash scenario of this paper be applicable? 

A special class of chondrites and chondrules is the CB 
and CH class. It has been recognized for some time that 
these chondrules may indeed have originated from collid¬ 
ing planetesimals or planets (see Desch et al. 2012 for a 
discussion). The evidence for this includes (1) the cryp¬ 
tocrystalline structure of CB chondrules, (2) their strong 
depletion of volatile elements, (3) zoned metal spherules 
that appear to have been condensed out of the gas phase 
in a highly energetic single-staged process (Krot et al. 
2005) and (4) their young age, 5 Myr after CAIs. One 
can interpret these properties as meaning (1) that these 
chondrules must have cooled more rapidly (seconds to 
minutes) than “normal” chondrules (hours) to explain 
their textures, (2) that the vapor-melt plume of the im¬ 
pact must have rapidly become of low enough density 
to allow volatile elements to get out of the chondrules 
and not recondense, (3) the impact must have been of 
very high velocity (several km/s) to produce metal vapor, 
(4) which appears consistent with the young age which 
means low gas densities in the protoplanetary disk and 
thus allows for high-speed impacts. Arguments (1) and 
(2) seem to require a rapid reduction of the density and 
optical depth, which can be achieved by a high impact 
velocity (also required for (3) and consistent with (4)) 
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and a low mass of the impactor. Krot et al. (2005) ar¬ 
gue, instead, for a very high mass collision to explain CB 
chondrules. Perhaps the adiabatic cooling of the impact- 
shock-compressed rock is enough to explain the fast cool¬ 
ing, or only a fraction of the mass from the giant impact 
actually gets converted into a plume of chondrules. 

In constrast to the case for the CB/CH chondrules 
our model may explain the properties of “normal” chon¬ 
drules: (1) barred or porphyritic textures indicating 
slower cooli ng (hours), (2) little depletion of volatiles 
(see Section [577]), (3) possibly lower temperatures and (4) 
earlier times after CAIs. The slower cooling rates come 
naturally out of our model because the optical depth of 
the cloud keeps the chondrules warm for a while. The 
little depletion of volatiles may be achieved by the high 
densities (though a quantitative analysis has to wait until 
paper II in this series). And the lower temperatures are 
natural for lower impact velocities which are expected in 
the earlier phases of the protoplanetary disk. The melt 
is the molten interior of the 26 Al-heated planetesimals, 
which also is expected during the earlier phases. 

We might also be able to apply our radiatively cool¬ 
ing expanding cloud model to the high impact velocity 
scenario for CB/CH chondrules, but it would require us 
to include more physics, in particular a proper equation 
of state for the shock-heated rock and the subsequent 
adiabatic expansion and adiabatic cooling. 

Among the “normal” chondrules there is also the issue 
of radial and barred textures versus porphyritic textures. 
The latter are the majority (about 85% of chondrules, 
Gooding & Keil 1981). Radial textures require that no 
nucleation sites were present in the melt, so that crys¬ 
tallization in the supercooled melt droplet starts from 
a single site and progresses radially outward from that 
site. In the splash scenario this means that the melt in 
the original pre-molten planetesimals was well above the 
liquidus, so that any previously existing crystals were 
molten. Conversely, for porphyritic chondrules numer¬ 
ous nucleation sites must have been present in the melt, 
and thus must have presumably remained present in the 
magma of the colliding planetesimals, meaning that that 
magma was likely slightly below the liquidus tempera¬ 
ture. The impact may have heated the melt through 
shock-heating above the liquidus, but if the cooling sets 
in quickly enough some of these nucleation sites may sur¬ 
vive. The question of porphyritic textures and the origin 
of the nucleation sites in the splash scenario remains, 
however, a tricky one. Sanders & Scott (2012) propose 
instead that these nucleation sites may have entered into 
the melt droplets as pollution from e.g. the surface re- 
golith that may have been mixed in with the chondrules. 
The problem with this (as with the relict grain issue of 
Section 15.511 is that the timing of mixing must be ideal: 
too early and the temperature is still above the liquidus 
and these regolith particles will also melt; too late and 
the chondrules have already formed radial or barred tex¬ 
tures. 

Finally there is the issue of Fe-content of chondrules: 


the H, L and LL classes. If the pre-molten planetesimals 
have experienced strong differentiation, it is possible that 
the mantles of these planetesimals are poor in iron and 
nickel. If the masses of the planetesimals are low enough 
and convection might be present, then perhaps all too 
strong loss of Fe in the mantle can be avoided. The issue 
of H, L and LL chondrites in the context of the splash 
scenario is further discussed in Sanders & Scott (2012). 


6 . SUMMARY AND CONCLUSION 

This paper studies the radiative cooling of an ex¬ 
panding cloud of lava droplets originating from the low- 
velocity collision between two partly or fully pre-molten 
planetesimals. The model is also valid for high-speed col¬ 
lisions between cold planetesimals, in which the melt is 
produced by the collision itself. 

The main result of this paper is a self-consistent cool¬ 
ing profile for chondrules based on the initial conditions 
of the expanding plume of melt droplets (chondrules). 
These cooling profiles are characterized by a time scale 
Aooi given by Eq. m and are well approximated by 
a constant temperature for t < t coo i and a powerlaw 
dropoff (Eq. [TH1) for t > t c 00 i- The cooling rate is given 
by Eq. (H3- The cooling time scale t coo \ depends on 
the total mass of the plume and on its expansion veloc¬ 
ity. Both parameters are related to the conditions of the 
collision between the two planetesimals. 

A large total mass of the expanding cloud of chondrules 
means that the colliding planetesimals must have been 
at least of that mass, but presumably considerably more 
massive since not all of the matter is likely to end up in 
a single expanding plume. The expansion velocity of the 
droplet cloud is related to the collision velocity of the 
two planetesimals: it is unlikely that it is higher than 
the impact velocity. Note that the expansion velocity 
is different from (and presumably much less than) the 
velocity by which the cloud moves away from the impact 
site. To relate both parameters (M c i ouc i and u e xp) to 
paramters of the colliding planetesimals it is necessary 
to perform 3-D hydrodynamic simulations of the impact, 
such as those performed by Asphaug et al. (2011). 

The results of our model, when compared to known 
cooling time scales of chondrules, put constraints on 
A/cioud and v exp , and thus on the conditions under which 
planetesimal collisions can produce chondrules. 
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APPENDIX 


TIME-DEPENDENT RADIATIVE TRANSFER: DISCRETIZATION AND IMPLICIT INTEGRATION 

Solving the radiative transfer moment equations coupled to the time-dependent heating/cooling equation requires 
an implicit integration scheme. The two coupled equations to solve are: 


Fi(J,T) := 


11 d (r 2 \d{fJ ) 3/ — 1 


J\> ~ J 


i(J,T) := pc v ^-j- — 47 tq - (j— —T A \ 

at V 7 r / 


The inner boundary condition is: 

Fi(J, T,r = r in ) ~ H(r = r in ) = 0 

where r; n is the smallest radius of our grid, which we take r; n <C i? c ioud- This condition translates into 

[ d{fJ) 3/ - 1 1 

+ -± - J = 0 

dr r J r=r . n 

The outer boundary condition is 

Fi(J, T,r = r 0 ut) := H(r = r ou t) - hJ = 0 

with h = -fffte(Tout)/<^fte(Tout), which translates into 

z WJ) , 3/ — 1 , ,1 


Fi(J, T, r = r 0 ut) := 


+ — -J + haJ 

dr r 


The only time-derivative in the equations is the one on the temperature. The radiation field is assumed to immediately 
adapt. 

We put J(r,t) and T(r,t) on a spatial grid , rjv}. So for time step n we have the values </”,••• , J 7 ^ 

and T", • • • ,Tjy. The above equations are also to be evaluated at these grid points: F\ \ = 0, • • • ,-Fi,jv = 0 and 
Ty i = 0, • • • , F 2 ,n = 0. Let us define 

Q n / jn rpn jn rpn jn rpn\T (\1\ 

~ \ J 1 ? 1 1 j J 2 i 1 2 i * * ' ? J Ni 1 N) J 


F(Q) = (F hl (Q),F 2 , 1 (Q),F h2 (Q),F 2 , 2 (Q),--- ,F hN (Q),F 2 , N (Q)Y 
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The objective of the time-integration of these equations is to find the values of Q n+1 : the values of Q at the next time 
step. 

For numerical stability we express all instances of Q in the equations F in their future form: Q" +1 , except when 
a time-derivative of Q is used, which is written as dQ/dt = ( Q n+l — Q n )j At and where the present value of Q n is 
required. We thus get as our set of equations: 

F(Q n+1 ) = 0 (A9) 

which we need to solve for Q n+1 . If F were linear in Q n+1 then one could write the vector F as a matrix multiplication 
with the vector Q n+1 plus a contant vector R: 


F = M Q n+1 + R 

(A10) 

Then the solution to F = 0 would be just a matrix inversion: 


Q n+1 = - M inv R 

(All) 

However, F is non-linear in Q n+1 because the Planck function B(T) = (a/ir)(T n+1 ) 4 . 
function F(Q n+1 ). So let us start with an initial guess . We typically then have 

And so we have a non-linear 

HQl ot 1 )^ 0 

(A12) 

Ideally we want to find a Q^t 1 for which F(QV’A 1 ) = 0, but we can only use Newton’s method by approximating 

F(Q"'+ ) using first order Taylor expansion: 


r)F 

F(Ql rt 1 ) - F (<3"o) ) + QQ ■ W& 1 - <?(o) l ) = 0 

(A13) 

or more in general for the fc-th iteration: 


F(Q^ 1)) - + ag • (Q( fe +1) ^ Qw 1 ) = 0 

(A14) 

If we define 

An»+i ,nn+i d™+ 1 

— b/( fc +i) 

(A15) 

we can solve for : 


+1 _ + i 

^ y nQ(k) ’ 

(A16) 

and then we obtain the new: 

n n+1 — AO n+1 + O n+1 

b/(fc+l) — 

(A17) 


Now let us write out the discrete equations Fi(J,T) and F 2 (J,T) explicitly: 


T n+i ) = 


1 


r i+ 1/2 fi+lJi+l — fiJi 


airf An \ a i+i/2 An+i/2 

r i- 1/2 fiJi - fi-iJi-i 


a i-l / 2 Al'i_ 1/2 


( 3 /i+1/2 — l)(Ji+l + Ji) r i+ 1/2 


2a. 


1 + 1/2 


(3/i- 1/2 — 1 1/2 


^ a i-l/2 


— J, d—(T" +1 ) 4 = 0 

7T 


For the gas equation we get: 


'yn+1 _ rpn 

F 2 ,i(J,T n+1 ) = pc v — 


At 


4710, (ji - ^(7j n+1 ) 4 ) — <? = 0 


Boundary conditions only have to be applied to the first equation. The inner boundary condition: 

F hl (J,T n+1 ) = /zJ2 ~ /lJl + J 2 ) 3f3/2 ~ 1 = 0 

a 3/2^ ? ’3/2 2 <+3/2 r 3/2 


(A18) 


(A19) 


(A20) 
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The outer boundary condition: 

Fi mj,t"+'> = + *-i) ( 3/y ~ 1/a ~ 1 

OiN-l/2^- r N-l/2 4 \ajV-l/2^JV-l/2 

The energy equation can be rescaled to: 

F 2 j(J,T” +1 ) = 7f +1 - If - (Ji - -(Tf +1 ) 4 ) - ^ = 0 

/9C V V 7T / /9C V 


+ h ) — 0 


The matrix coefficients then become: 

( 9Fi,i \ 1 , | 

dJi ) otivfAri t a i _i/ 2 Ar i _ 1 /2 a. i+1 / 2 Ar i+1 / 2 

1 


r 2 
r i- 1/2' 


:/< 


r 1+1/2^ 


- 1 


/ \ _ 1 f r i+l/2/ i + 1 (3/i+l/2 — l)?'i+l/2 1 

V9J i+ i/ aurfAri | a i+ i/ 2 Ar i+ i/2 2a i+1 / 2 j 

/ cfFqj \ _ 1 f r i-l/2fi~ 1 (3/i—1/2 ~ l)r*-l/2 1 

\dJi-i) our?Ari y a i _i/ 2 Ar,_i/2 2a l _ 1/2 J 

/ c>F M \ 4a +1 3 

V577+ 1 / Jr 1 * j 

/ _ \-KUj 


\ d Ji J PiCy 

-K" + ‘) 3 A* 
' v aT/ i+ v 7T 

The matrix coefficients for the boundary conditions become: 

/dFi.A 1 f 3/3/2-1 | 

\ 9 J 2 J 0:3/2 l 2r3/2 

/ dF 1 , 1 \ 1 f 3 / 3 / 9 -1 

I r\ 


/2 

Ar 3/2 

h 


J _ ■/! 

O 3/2 1 2r3/2 Ar 3 / 2 

1 / 3/at_i/2 — 1 


v 9Ji; 

( dF N>1 \ = _ , , 

\ 9Jat / ojv_i/2 V 2rjv~i/2 
/ aFjV.l V. 1 / 3/AT-1/2-1 _ 

\8Jn-iJ Ojv-i/2 V 2rjv~i/2 Arjv-1/2 


/w 

A r N _i/ 2/ 
/jv-i 


ft 

2 

ft. 

2 


(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 

(A31) 

(A32) 


OjV —1/2 



